航天器视场遮挡分析方法研究

敏感器视场遮挡分析是航天器构型布局设计的重要组成部分，针对大型在轨变构型航天器敏感器视场遮挡分析工作量大、出错率高、遮挡率无法量化等特点，本文提出一种基于敏感器第一视角的视场遮挡快速分析方法。该方法通过设置第一视角图像蒙版及图像映射，模拟方锥形、圆锥形、球形等各类敏感器视场，通过逐个像素比对，计算视场遮挡率和遮挡方位。使用此方法建立了空间站敏感器视场遮挡分析系统，结果表明新方法能够显著提高视场遮挡分析效率和准确率，并且为机构类部件运动过程对敏感器视场遮挡的分析提供了有效工具。     

非线性能量阱对飞轮振动抑制效果的实验研究

利用薄片金属梁作为非线性刚度元件，构建一种非线性能量阱。为了验证非线性能量阱对飞轮振动特性的抑制效果，搭建地面飞轮测试系统。分别采用力锤敲击和飞轮以恒定转速运行的方法施加不同载荷，对安装非线性能量阱前后测试系统各典型位置的振动响应进行对比研究。利用实测数据，采用傅里叶变换方法将动态响应信息分别在时间域和频率域内展开。研究表明：力锤敲击产生的冲击载荷和飞轮以恒定转速运行产生的稳态载荷均可以激发非线性能量阱定向能量传输；受非线性能量阱的影响，柔性支架设备安装面和飞轮安装面振动响应得到有效抑制。但是，由于非线性能量阱阻尼和吸振子质量较小，不同程度削弱了其振动抑制效果。     

飞机液压系统柱塞泵热分析

随着飞机液压系统向高压化、大功率化方向发展，液压系统的发热及温升问题受到了广泛关注。柱塞泵作为液压系统中重要的能量转换装置，对液压系统的功率损失和温度有重大的影响，良好的柱塞泵的热分析对液压系统的热设计和温度控制有重要的意义。然而传统的柱塞泵热分析方法，如平均油温法、软件仿真法和热节点网络法，因不够准确、高效和简洁而不能满足现代飞机液压系统热分析和热设计的相关需求。本文针对飞机液压系统柱塞泵发热的相关问题，采用基于集总参数的热网络分析方法，建立了柱塞泵热网络模型，通过对柱塞泵不同工况下的仿真，验证了模型的有效性和准确性，为柱塞泵及液压系统的热分析提供了参考。     

空间载荷机柜热控设计的仿真研究

空间站主动热控系统的作用是提供和保持航天员、仪器设备和结构部件所要求的温度、湿度和风速等热环境。本文以空间载荷机柜为热控系统的基本单元，建立了以单相流体循环回路为基础、采用液冷结合风冷的冷却方式的载荷机柜数学模型；介绍了模糊增量控制算法的工作原理；提出了采用模糊增量控制算法，以液冷支路进出口温差和风冷支路冷却量占总冷却量的比值为控制目标的反馈控制策略。同时文章采用数值模拟方法研究了控制系统的动态响应效果。仿真结果显示，模糊增量控制算法具有响应迅速、超调量小、无静态误差的特点，控制策略满足有效性、稳定性和高精确度的要求。     

NiCrAl/Diatomite可磨耗封严涂层制备及性能

分别利用大气等离子喷涂技术制备了Ni/Al粘结底层、火焰喷涂技术制备了NiCrAl/Diatomite可磨耗封严涂层，研究了喷涂距离、送粉速率、火焰气体总流量、喷涂角度和氧燃比等喷涂参数对生长速率和涂层硬度的影响。结果表明，在较低的氧燃比条件下（O₂∶C₂H₂<1.6），随着喷涂距离的增大，涂层生长速率逐渐下降，硬度先保持不变后逐渐上升；当O₂∶C₂H₂>1.6时，涂层的生长速率会随着喷涂距离的增大先增大后减小，硬度随之逐渐下降。随着喷枪对基体相对移动速度的增大，涂层生长速率略微下降，硬度略有上升。气体总流量的增大使得涂层的生长速率明显上升，硬度明显下降。随着送粉量的增大，涂层生长速率明显提高，硬度随之先上升后下降。     

Comparison of GPS derived TEC with the TEC predicted by IRI 2012 model in the southern Equatorial Ionization Anomaly crest within the Eastern Africa region

We have compared the TEC obtained from the IRI-2012 model with the GPS derived TEC data recorded within southern crest of the EIA in the Eastern Africa region using the monthly means of the 5 international quiet days for equinoxes and solstices months for the period of 2012 – 2013. GPS-derived TEC data have been obtained from the Africa array and IGS network of ground based dual-frequency GPS receivers from four stations (Kigali (1.95°S, 30.09°E; Geom. Lat. 11.63°S), Malindi (2.99°S, 40.19°E; Geom. Lat. 12.42°S), Mbarara (0.60°S, 30.74°E; Geom. Lat. 10.22°S) and Nairobi (1.22°S, 36.89°E; Geom. Lat. 10.69°S)) located within the EIA crest in this region. All the three options for topside Ne of IRI-2012 model and ABT-2009 for bottomside thickness have been used to compute the IRI TEC. Also URSI coefficients were considered in this study. These results are compared with the TEC estimated from GPS measurements. Correlation Coefficients between the two sets of data, the Root-Mean Square Errors (RMSE) of the IRI-TEC from the GPS-TEC, and the percentage RMSE of the IRI-TEC from the GPS-TEC have been computed. Our general results show that IRI-2012 model with all three options overestimates the GPS-TEC for all seasons and at all stations, and IRI-2001 overestimates GPS-TEC more compared with other options. IRI-Neq and IRI-01-corr are closely matching in most of the time. The observation also shows that, GPS TEC are underestimated by TEC from IRI model during noon hours, especially during equinoctial months. Further, GPS-TEC values and IRI-TEC values using all the three topside Ne options show very good correlation (above 0.8). On the other hand, the TEC using IRI-Neq and IRI-01- corr had smaller deviations from the GPS-TEC compared to the IRI-2001.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Simulating planetary wave propagation to the upper atmosphere during stratospheric warming events at different mountain wave scenarios

Parameterization schemes of atmospheric normal modes (NMs) and orographic gravity waves (OGWs) have been implemented into the mechanistic Middle and Upper Atmosphere Model (MUAM) simulating atmospheric general circulation. Based on the 12-members ensemble of runs with the MUAM, a composite of the stratospheric warming (SW) has been constructed using the UK Met Office data as the lower boundary conditions. The simulation results show that OGW amplitudes increase at altitudes above 30 km in the Northern Hemisphere after the SW event. At altitudes of about 50 km, OGWs have largest amplitudes over North American and European mountain systems before and during the composite SW, and over Himalayas after the SW. Simulations demonstrate substantial (up to 50–70%) variations of amplitudes of stationary planetary waves (PWs) during and after the SW in the mesosphere-lower thermosphere of the Northern Hemisphere. Westward travelling NMs have amplitude maxima not only in the Northern, but also in the Southern Hemisphere, where these modes have waveguides in the middle and upper atmosphere. Simulated variations of PW and NM amplitudes correspond to changes in the mean zonal wind, EP-fluxes and wave refractive index at different phases of the composite SW events. Inclusion of the parameterization of OGW effects leads to decreases in amplitudes (up to 15%) of almost all SPWs before and after the SW event and their increase (up to 40–60%) after the SW in the stratosphere and mesosphere at middle and high northern latitudes. It is suggested that observed changes in NM amplitudes in the Southern Hemisphere during SW could be caused by divergence of increased southward EP-flux. This EP-flux increases due to OGW drag before SW and extends into the Southern Hemisphere.     

Comparison of vertical E × B drift velocities and ground-based magnetometer observations of DELTA H in the low latitude under geomagnetically disturbed conditions

In the present work, we analyzed the daytime vertical E × B drift velocities obtained from Jicamarca Unattended Long-term Ionosphere Atmosphere (JULIA) radar and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca, and Piura in Peru for 22 geomagnetically disturbed events in which either SC has occurred or Dst_max < ?50 nT during the period 2006–2011. The ΔH component of geomagnetic field is measured as the differences in the magnitudes of horizontal H component between magnetometer placed directly on the magnetic equator and one displaced 6–9° away. It will provide a direct measure of the daytime electrojet current, due to the eastward electric field. This will in turn gives the magnitude of vertical E × B drift velocity in the F region. A positive correlation exists between peak values of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of the events. It was observed that 45% of the events have daytime vertical E × B drift velocity peak in the magnitude range 10–20 m/s and 20–30 m/s and 20% have peak ΔH in the magnitude range 50–60 nT and 80–90 nT. It was observed that the time of occurrence of the peak value of both the vertical E × B drift velocity and the ΔH have a maximum (40%) probability in the same time range 11:00–13:00 LT. We also investigated the correlation between E × B drift velocity and D_st index and the correlation between delta H and D_st index. A strong positive correlation is found between E × B drift and D_st index as well as between delta H and D_st Index. Three different techniques of data analysis – linear, polynomial (order 2), and polynomial (order 3) regression analysis were considered. The regression parameters in all the three cases were calculated using the Least Square Method (LSM), using the daytime vertical E × B drift velocity and ΔH. A formula was developed which indicates the relationship between daytime vertical E × B drift velocity and ΔH, for the disturbed periods. The E × B drift velocity was then evaluated using the formulae thus found for the three regression analysis and validated for the ‘disturbed periods’ of 3 selected events. The E × B drift velocities estimated by the three regression analysis have a fairly good agreement with JULIA radar observed values under different seasons and solar activity conditions. Root Mean Square (RMS) errors calculated for each case suggest that polynomial (order 3) regression analysis provides a better agreement with the observations from among the three.     

Role of Equatorial Anomaly in Earthquake time precursive features: A few strong events over West Pacific zone

The earthquake (EQ) time coupling processes between equator-low-mid latitude ionosphere are complex due to inherent dynamical status of each latitudinal zone and qualified geomagnetic roles working in the system. In an attempt to identify such process, the paper presents temporal and latitudinal variations of ionization density (foF₂) covering 45°N to 35°S, during a number of earthquake events (M?>?5.5). The approaches adopted for extraction of features by the earthquake induced preparatory processes are discussed in the paper through identification of parameters like the ‘EQ time modification in density gradient’ defined by δ?=?(foF_{2 max}???foF_{2 min})∕τ_mm, where τ_mm – time span (in days) between EQ modified density maximum and minimum, and the Earthquake time Equatorial Anomaly, i.e. EEA, one of the most significant phenomenon which develops even during night time irrespective of epicenter position. Based on the observations, the paper presents the seismic time coupling dynamics through anomaly like manifestations between equator, low and mid latitude ionosphere bringing in the global Total Electron Content (TEC) features as supporting indices.